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Abstract 



The photo-production of a pair of scalar particles in the presence of an intense, circularly polarized laser beam is investigated. 
Using the optical theorem, explicit expressions are given for the pair production probability in terms of the imaginary part of 
the vacuum polarization tensor. Its leading asymptotic behavior is determined for various limits of interest. The influence of the 
absence of internal spin degrees of freedom is analyzed via a comparison with the corresponding probabilities for production of 
C\| .spin-72 particles; the lack of spin is shown to suppress the pair creation rate, as compared to the predictions from Dirac theory. 
(^Potential applications of our results for the search of minicharged particles are indicated. 
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1. Introduction 

q_. Understanding the nonlinear and unstable nature of the quan- 
1 turn vacuum in the presence of a strong electromagnetic field 
^^constitutes an important task of theoretical physics. Correspond- 
(— 1 ^ing studies have revealed a nontrivial vacuum structure, suit- 
1— i-able to explore the low-energy frontier of particle physics 0], 
^ (111- Moreover, perspectives of achieving ultrahigh field intensi- 
^ -ties (/ ~ 10 26 W/cm 2 ) in short laser pulses of few femtosec- 
ly-j [onds duration (J 01 have motivated a growing interest in the 
Q\ phenomenology purely associated with the quantum nature of 
l/~) the electromagnetic interaction (see JHIfJ for recent reviews). 
This is because the envisaged laser field strengths lie only 1- 
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2 orders of magnitude below the critical value E c = 1.3 x 
10 16 V/cm where QED vacuum nonlinearities become substan- 
tial and spontaneous vacuum decay into electron-positron (e~e + ) 
airs via the Schwinger mechanism is expected to occur 



In combination with an incident high-energy particle, strong 
laser fields can induce e~e + pair production already at inten- 
sities available today. In a pioneering experiment at SLAC 
Jioll . a multi-GeV photon decayed into a pair while propagat- 
ing through a moderately intense laser pulse (/ ~ 10 18 W/cm 2 ). 
This process, involving the simultaneous absorption of several 
laser photons, represents a nonlinear version of the well-known 



Breit-Wheeler reaction 111 Illl2lll3ll . The high-energy non-laser 



photon originated from Compton backscattering of SLAC's ul- 
trarelativistic electron beam off the laser pulse. In the near fu- 
ture, corresponing studies can be conducted within all-optical 
setups using laser-accelerated relativistic electrons as projec- 
tiles [14]. Other pair production mechanisms may be probed in 
ultrarelativistic proton-laser collisions 11151 llq, Il7l Il8l Il9ll . 



In view of the upcoming high-field laboratories UQ], theo- 
reticians are currently investigating further properties and appli- 
cations of photo-induced e~e + pair production in intense laser 
fields. For example, due to their broad frequency composition, 
laser pulses of ultrashort duration have been shown to modify 
the created particle spectra lEoll and lead to characteristic en- 
hancements in the pair production probability 12 ll 12211 . Photo- 
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induced pair production also plays a crucial role for the devel- 
opment of QED cascades which may give rise to e~e + plasmas 
of very high density J23 , 24 ] . The nonlinear Breit-Wheeler pro- 
cess moreover offers a promising means to measure ultrashort 
y-ray pulses via e~e + streaking 12511 . Superimposing the field of 
a high-energy photon onto a strong electric field may also help 
catalysing the Schwinger effect 

In the present Letter, we study the photo-induced creation of 
a pair of spin-0 particles in the presence of a strong monochro- 
matic laser beam. Our motivation is twofold. First, while the 
probabilities for the creation of fermion pairs are known for a 
long time Il2lll3ll . it is relevant to establish the corresponding 
formulas for scalar particles because they can be useful for the 
ongoing search of minicharged particles which may have ei- 
ther fermionic or bosonic character iB Esll . Second, our results 
provide insights into the fundamental question as to how the 
spin degree of freedom affects the photo-induced pair produc- 
tion process. To this end, a comparison with the known results 
for fermion pair production will be drawn. Such an informa- 
tion complements previous works where spin-resolved calcula- 
tions of the nonlinear Breit-Wheeler process via helicity ampli- 
tudes [290 an d the internal spin polarization vector OOll were 
performed. We note besides that comparative studies between 
the behaviour of bosonic and fermionic particles in strong laser 
fields have rececently been carried out with respect to Comp- 
ton, Mott and Kapitza-Dirac scattering J5, 31, 32], nonlinear 
Bethe-Heitler pair creation in proton-laser collisions [33] and 
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the Klein paradox (34 , 351. 

Our theoretical approach relies on the polarization tensor 
Il /JV (A:i , k 2 ) of scalar particles in the presence of a strong laser 



field M 13711 whose imaginary part is related to the pair pro- 
duction probability via the optical theorem. While the polariza- 
tion tensor for Dirac fermions has already been exploited suc- 
cessfully to calculate various e~e + pair production processes in 
strong laser fields 1 171 , 1 181 13611 . to the best of our knowledge the 
present calculations represent the first application of the corre- 
sponding polarization tensor for the scalar case. 

2. General considerations 

To begin with, let us consider the field of a plane electro- 
magnetic wave of the forrrQ 



x) — a^Vi(xrx) + c^i 2 (xx), 



(1) 



with Ai 2 denoting the wave amplitudes and 1/0,2 being arbitrary 
functions. The wave four-vector x? = {x°, x) fulfills the rela- 
tions x 2 = and xa\ = xa 2 — a\a 2 — 0. According to l36l 371. 
the vacuum polarization tensor in this field, 

= ci Kfhf + c 2 hf A« + csA^Af' 
+ c 4 A^At 2) +C5 A«At 4) (2) 

can be expanded in terms of a basis set of Lorentz covariant 
vectors A® which are constructed from fundamental symmetry 
principles. They are explicitly given by 
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- x^cCy' - x v aj^ (i = 1,2) are the amplitudes 
of the external field modes. We note that the short-hand no- 
tation k in Eq. 0]) may stand for either k\ or k 2 . It is worth 
mentioning at this point that, for k — k\, the vectors Ai, A2 
and A3 are orthogonal to each other, A^ A^ = —S l i, and ful- 
fill the completeness relation g^ - -p^ 1 - - A® A® with 
gfjy = diag(+l, -1,-1,-1) denoting the metric tensor. A sim- 
ilar statement applies if the set of vectors Aj, A2 and A4 with 
k — k 2 is considered. We emphasize that Eq. (0 does not de- 
pend on which choice of k is taken since the difference between 
£1 and k 2 is proportional to x [see Eq. (0 below]. 

The form factors c, in Eq. (f2]i are distribution-valued func- 
tions which depend on the field shape via the functions fr. They 
have been evaluated thoroughly for the case of spin-^ particles 
in 1 3(1 37]. Also for the case when the virtual charge carriers 
in the Feynman loop are spin-0 particles general expressions 
for the Cj were provided in these references; but these formulas 
were not further evaluated. 



Using the general expressions from 1361 13711 and assuming 
that the laser field is an elliptically polarized wave with 



iffi = cos(xx) and \fi 2 = sin(>a:), 



(4) 



we find that the form factors in Eq. © for the scalar case are 
given by 
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Here, a = e 2 is the fine structure constant, e and m denote the 
particle charge and mass, respectively, and 



xk 



* = ^Z?' 6 = ~ 



a = 1,2). 
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As Eq. (|5j shows, the polarization tensor decomposes into elas- 
tic (k\ = k 2 ) and inelastic {k\ + k 2 ) parts. Those terms which 
contain the Dirac deltas 5 A {k\ — k 2 + 2Nx) with N + are 
responsible for the inelastic scattering of a photon in the field 
of the wave. For our purposes, however, only the elastic part 
is relevant. The corresponding functions d^°\ i = 1,2,3,4,5 
contained in Eq. (f5]l are given by 

df = -df = ^ 2 pA J Q (z)sign[A], 

1 f 2 
df = --f x A, (j (z) - iJ' (zj) + f sm 2 (p)J (z) 



i\A\ (l - v 2 ) 



8p 
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where Jo(z) is the Bessel function of zero order and J' (z) its 
derivative. The remaining parameters are 
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and A\ = A + 2A . 

A substantial simplification is achieved when the external 
field is taken as a circularly polarized wave (£1 = % 2 = £). In 
this case, we find it convenient to express the elastic contribu- 
tion as 



n^Ofci.fc) = /(27r) 4 5 4 (^! - k 2 )n„ v (k 2 ) 
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with 



IVfe) = (A<X 2 > - A®A?>)af + (A^Af + A® A®) 

(10) 



'From now on "natural" and Gaussian units c = h = Ane® = 1 are used. 
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Figure 1: Diagrammatic representation of the optical theorem applied to the 
photo-production of a pair of spinless particles in the field of a wave. In the 
left-hand side we have represented the squared modulus of the the S -matrix 

element integrated over the final phase volume f, = f ^-^t . Here the 

° r J J J (2ny (2ny 

doubles lines refer to the exact Klein-Gordon states interacting with the exter- 
nal background. The right-hand side, however, contains the imaginary part of 
the vacuum polarization tensor where the double lines represent the exact prop- 
agators in the field of the wave. In both sides the wavy lines denote the legs 
corresponding to a photon field. 



Here the involved coefficients are given by 
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Of> = 2£ 2 A psign [A] , nf = -^v 2 (1 - e>) , 
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Now, the unitarity condition of the dispersion S -matrix pro- 
vides the optical theorem. According to the latter, the total 
creation rate of a pair of spin-0 particles from a real photon 
(k 2 = 0) with polarization vector turns out to be 9^ 0) = 
e£e* f Im T¥' v /to where a> is the photon frequency (see Fig. [TJ. 
If the photon is unpolarized, one has 
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where the completeness relation g MV = — Yit ^i e T was use( l m 
the second step. Because of this fact, the expression above re- 
duces to 



= Imnf 



So, among the set of form factors previously described in Eqs. 
(fTTI) - (fT2l) . only one contributes to CH (0) . In correspondence, the 
probability of producing a pair of scalar particles by a photon 
travelling through a laser field of circular polarization can be 
represented by the following double-parametric integral 

2ttw Jo Jo p 1 
+ 2£ 2 sin 2 (p)}, (15) 

where the abbreviation 77 = {pi A) [l + f 2 (l - sin 2 (p)/p 2 )] has 
been used. Note that for a real non-laser photon, the parameter 
A in Eq. ((6]) is always positive. 

The structure of Eq. (fT~5T > shares certain similarities with 
the corresponding production rate of spin-^ particles (see 



Eq. (D.19) in 0711 ). It differs, however, by (i) an overall factor 
of 2 which coincides with the spin-— multiplicity, and (ii) the 
precise form of the integral over v. As we will see, this gives 
rise to differences between 5H (0) and SH's' regarding both the ab- 
solute size and the functional dependence. We should mention 
at this point that an exact evaluation of 9^ (0) is quite difficult 
to perform. However, the dependence on the parameters of the 
theory, i.e. A and allows to obtain closed-form analytical rate 
expressions in various asymptotic regimes of interest. They will 
be derived in the forthcoming section. 



3. Asymptotic regimes and comparison with spin- 1 case 

Before undertaking the main calculations of this section, we 
briefly analyze the kinematics associated with the decay pro- 
cess y{k) + njiiyc) — > f + £ + , where n denotes the number 
of participating laser photons ji. In the center-of-mass frame, 
the corresponding energy-momentum balance implies the con- 
dition nkx = 2s 2 , with s being the (laser-dressed) energy of a 
final particle state. Accordingly, whenever the number of ab- 
sorbed laser photons exceeds the threshold value 

2m 2 t 

n* = - — (16) 
kh 

the decay can occur. Here, m„ = m(l + £ 2 )^ 2 denotes the ef- 
fective mass of the particle. It comes out as a consequence 
of considering the dressed four-momentum q M - p M + j^m^ 
(with p M the free four-momentum of the particle, i.e. p 2 = m 2 ) 
as the kinematical variable involved in the energy-momentum 
conservation jl2lll3ll . We point out that in the center-of-mass 
frame, the relative speed between the created particles is given 
by |#reil = |f- - f+l = 2v with 



P = = (l - — f 



(17) 



3.1. Two photon-reaction at £, < 1 

We wish to specialize Eq. (fT3T > to the case where £ < 1. 
In this context, we note that the probability of creating a pair 
is suppressed if the condition A «: 1 holds, so that dispersive 
processes become of main interest her^|. On the contrary, if 

(14) 

A ^ 1 + £ 2 (i.e. n* ^ 1), the pair production could take place 
with the absorption of just one photon from the laser wave (two- 
photon reaction). In this case, which we shall consider now, the 
oscillatory contribution present in the exponent of (fT~5T > becomes 
very small in comparison with the remaining terms, allowing to 
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variable according to p 
(I151 l. we arrive at 



without obstruction. After a change of 
p( 1 + £ 2 ) in the second term of Eq. 
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2 In the aforementioned limit, the theoretical description of such phenomena 
can be performed by means of an effective Lagrangian approach. For details 
we refer the reader to (^.[39|,[40|.[4l|l. 
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In order to provide an explicit expression of £H (0) we inte- 
grate by parts the terms containing the factor proportional to 
1 /p 2 . Afterwards, the residue theorem is applied to obtain 
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where 6(x) denotes the unit step function. The remaining inte- 
gral over v can be taken analytically without complications. It 
leads to 
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with «* given in (fT6l l. The respective expression for JHCi) can 
be read off from Eq. (Got by multiplying the latter by -2 and 
replacing the coefficient in front of the logarithmic function by 
1 + ««/(l + f) - «;/2(l + £ 2 ) (see Appendix D in Q). 

Let us consider some limiting cases. At «„ ^ 1, i.e. in the 
nonrelativistic limit (v <k 1), Eq. ( f20T > behaves as 
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whereas So, in this limit, the ratio between the 

production probabilities coincides with the spin multiplicity of 
a Dirac particle. Incidentally, this relation is also manifest to 
leading order between the rates of the respective Schwinger 
mechanisms 14211 . In contrast, for n, ~ 0, the created parti- 
cles are ultrarelativistic (v ~ 1) and the probability becomes 
independent of A : 
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The situation is quite different when the photo-production of 
Dirac particles is considered. In fact, under the same circum- 
stances, the creation rate of a pair of spin-i particles has a log- 



arithmic dependence: -=-£- [21n(2e/m«) — 1] ; see also 

1 3811 . As a consequence, the ratio between the fermionic and 
scalar rates grows logarithmically, i.e. D^W/OV ' ~ 41n(2s/m„) 
as s ~ (kx) 1 ^ 2 » m*. Such a behaviour resembles the cor- 
responding result arising in Compton scattering. Indeed, as is 
well known, when the energy of the incoming photon is very 
large (to » m) the total Compton cross section computed for 
Dirac fermions turns out to be (A J ) ~ — [ln(2w/m)+ 1/21. In 

com L v ' y ' J 

constrast, the leading asymptotic behaviour of the correspond- 
ing total cross section determined for spin-0 particles reads cr ' « 



2a 2 k 



3ill . Accordingly, cr^/c 



.to) 



■\n(2to/m) as to » m. 
Obviously, in the case under consideration, can be 

substantially larger than the statistical factor of 4 associated 
with the possible spin configurations in the final state. 

The results obtained in this section are summarized in Fig. [2] 
which displays the ratio DtCO/JH® as a function of A. Several 
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Figure 2: Relative dependence of the femionic and scalar pair production rates 
on the parameter A for different values of f < 1 , as indicated. The case f <8C 1 
corresponds to the Born approximation. 



curves are shown corresponding to different values of ^. The 
picture also includes the limiting case of f « 1 (solid line), 
which is compatible with the Born approximation. The explicit 
expressions associated with this limit can be read off from Eqs. 
d20b-d2Tb by replacing «, — > «b = l/A and setting £ — within 
the squared brackets. The latter procedure leaves us with a 
quadratic dependence on the parameter Moreover we note 
that, at larger values of ^ as 1 and A ~ 1 + £ 2 , next-to-leading or- 
der terms become increasingly important (see also lfl7lo . These 
terms which may give some minor contribution to the rates for 
£ = 1/2 have not been included in Fig. [2] 

3.2. Leading behaviour of '9^°' at asymptotically large % » 1 

In order to undertake the calculations in the high-field do- 
main with £ » 1 it is convenient to express the integral of the 
first two terms in Eq. ( fl3b in a more suitable form: 
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with Aq given in Eq. ([8]). Inserting this expression into Eq. ( fl31 ) 
we obtain 
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where the lower boundary of the p-integral was extended to -oo 
taking into account the symmetry of the integrand. Besides, the 

change of variable u = (l - v 2 ) has been carried out. We 
point out that the integral over this variable does not diverge, 
although the integrand is a singular function at u = I. Since the 
latter, moreover, decreases very fast like ~ 1/m 2 when u — > oo, 
it is expected that the main contribution results from the region 
around u ~ 1 . The situation is somewhat different with respect 
to the variable p. While in this case the integrand falls off for 
p — » +oo as well, it is a regular function in p with vanishing 
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limit at p — » 0. In order to elucidate the mainly contributing 
region of this integral, we split the integration domain into three 
regions: from -oo to -po (lower region), from -po to po (inner 
region) and from po to oo (upper region). This is, 



r° dp _ r Po dp r p ° dp r 

xJ — OO P xj — OO P xj — pQ P xj Pq 



dp 
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where the positive dimensionless parameter po is chosen to ful- 
fill the conditions 



£ «: po «: 1 and 



(A/f)' /3 «p Q « 1. 



(25) 



Within the inner narrow integration region, where \p\ ^ po «c 1, 
we can Taylor expand the respective integrand and, afterwards, 
perform a change of variable according to s — p£. The latter is 
also carried out in the integrals over the lower and upper regions 
but, after extending the resulting integration limit po£ — » oo, no 
relevant contribution comes from them. Thus, the total integral 
over p can be well approximated by 
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We remark that Eq. d26l is accurate up to terms that decrease 
exponentially, like ~ (po£) _1 exp [ _ x|=(po£) 3 ] or even faster. 

The integration over s can, then, be done with the help of 
the Macdonald function K v (x) using the relational 
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We insert the resulting expression into Eq. d24l and use the 
identity -xKi /3 (x) = x£Ky 3 (x) + jK2/ 3 (x) (see Eqs. (8.486.12) 
and (8.486.16) in JH) to obtain 
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where, in addition, an integration by parts has been carried out. 
Besides, the abbreviation £ = A£ — J-(l-£-£) has been intro- 
duced, with the critical field strength E c = m z /e. The parameter 
( encodes the two paths to vacuum polarization effects: either 
by increasing the external field amplitude or by increasing the 
photon energy. It also depends on the propagation directions 
of the photon k = k/\k\ and the strong laser field k = x/\x\- 
Observe that £ is maximized when the fields counterpropagate, 
i.e., whenfc ■ k = -1. 



^Equations (3.695.1-2) in Ref. allow to find the identity 
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After differentiating with respect to b and using Eqs. (8.486.10-1 1) of the afore- 
mentioned reference we establish the relations j27M28t. 



We recall that Eq. (1291 was derived under the assumptions 
given in d25l . Thus, the rate expression above applies when 
the number of absorbed laser photons is very large, n n* ~ 
% 2 /A » 1. In this case, many photon orders n contribute to 



the production rate 112111311 . Note that Eq. d29l is independent 
of the frequency of the strong brackground field. The prob- 
lem thus becomes quasi-static with respect to the external wave 
which, consequently, may be approximated by a pure constant 
crossed field. We have checked that, starting from the general 
expression of the vacuum polarization tensor Tl^ 1 * in such an 



external field configuration 14411 . we can also arrive at Eq. ( 1291 . 

The corresponding production rate of spin-^ particles emer- 
ges from ( |29l via the replacement: (4m - 1) — » 2(8m + 1). The 
resulting integrand turns out to be greater than the respective 
one of Eq. d29l within the whole integration domain of the 
M-variable. This fact guarantees that > 1. Hence 

the internal spin degrees of freedom promote the creation rate 
of spin-i particles as compared with scalar ones. 

So far, no restriction has been imposed on the parameter £. 
Let us consider the situation in which £ » 1 . To be consistent 
with our previous condition (£ » 1) the parameter A must be 
restricted to A » 1 /£. We may exploit the small-argument be- 
haviour of the Macdonald functions, K v (x) ~ ^ 0)" 0], and 
obtain 
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with Y(x) denoting the Gamma function. The corresponding 
result for the fermionic case reads 13711 



5am 2 ^ (3 
6 V^t, 



(31) 



Thus, the role of the internal spin degrees of freedom are man- 
ifest as 9t(i) =59\ (0) . 

On the contrary, if ( <k 1 (corresponding to A «: 1 at 
£ » 1), the large asymptotic behaviour of the Macdonald func- 
tions applies, i.e. K v (x) ~ yp^e~ x 114311 . With this expansion 
in mind, the integration over u can be performed. The latter 
becomes particularly simple because the region u ~ 1 provides 
the essential contribution. This leads to 
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A comparison with the corresponding fermionic rate 13811 
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(33) 



leads to write . It is interesting to note that a sim- 

ilar result was found for the strong-field Bethe-Heitler process 
ycon\+njL — > ■ In the parameter domain where £ » 1 and 

E ~ E c , the rate associated with the creation of Dirac fermions 
in the Coulomb field of a nucleus exceeds by a factor of ^ 7 the 
corresponding rate for spin-0 particles [33]. 
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4. Conclusions and outlook 

Photo-initiated production of a pair of spinless particles in 
the field of a circularly polarized laser wave was investigated. 
Compact expressions for the pair production rate were obtained 
in various asymptotic parameter regimes. Our analysis was car- 
ried out by considering the imaginary part of the vacuum polar- 
ization tensor as dictated by the optical theorem. Comparisons 
between the creation rates for spin-^ versus spin-0 particles re- 
vealed how spin effects are manifest at different energy scales. 
It was shown that the spin degrees of freedom affect the ab- 
solute magnitude and, in the limit of high photon energies, also 
the functional form of the production rates. The rate predictions 
based on the Dirac theory were generally found to be signifi- 
cantly larger than the corresponding results for Klein-Gordon 
particles. 

We note that the absence of internal spin degrees of free- 
dom in the scalar theory renders the evaluation of the vacuum 
polarization tensor somewhat easier than in the fermionic case. 
It is worth mentioning that, since the total production rates for 
spin-0 and spin-^ particles in the high-field limit » 1) dif- 
fer by an overall factor only, this technical simplification could 
be exploited to extract new insights into the production process 
which are valid not only for scalar particles but also for e~e + 
pairs. Indeed, the very similar behaviour of the production rates 
in this regime indicates that the process is mainly determined by 
the features of the weak photon and strong laser fields, rather 
than the particular property of the matter field. 

Finally we point out that, whenever the energy scale re- 
mains within the phenomenological limits of QED and its fun- 
damental principles are preserved, the photo-production of a 
hypothetical pair of spinless particles characterized by a tiny 
fraction e of the electron charge —\e\ would not differ qualita- 
tively from the respective creation of a spinless electron-positron 
pair. In correspondence, the rate associated with the latter phe- 
nomenon can be obtained from the expression derived in this 
Letter by replacing the electron parameters (|e|,m) by the re- 
spective quantities (e|e|, m c ) associated with a minicharged par- 
ticle. Therefore, our results can be useful in the ongoing search 
for minicharged particles |jl 0, 0] and help to improve our 
understanding of how the spin degrees of freedom affect the 
relevant experimental observables. 
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